arXiv: 1507.07161 v3 [cs.NI] 18 Oct 2015 


l 


An Optimal Resource Allocation with Frequency 
Reuse in Cellular Networks 

Ahmed Abdelhadi and T. Charles Clancy 
Hume Center, Virginia Tech, Arlington, VA, 22203, USA 
{aabdelhadi, tec}@vt.edu 


Abstract —In this paper, we introduce a novel approach for 
optimal resource allocation with frequency reuse for users 
with elastic and inelastic traffic in cellular networks. In our 
model, we represent users’ applications running on different user 
equipments (UE)s by logarithmic and sigmoid utility functions. 
We applied utility proportional fairness allocation policy, i.e. 
the resources are allocated among users with fairness in utility 
percentage of the application running on each mobile station. 
Our objective is to allocate the cellular system resources to 
mobile users optimally from a multi-cell network. In our model, 
a minimum quality-of-service (QoS) is guaranteed to every user 
subscribing for the mobile service with priority given to users 
with real-time applications. We show that the novel resource 
allocation optimization problem with frequency reuse is convex 
and therefore the optimal solution is tractable. We present a 
distributed algorithm to allocate the resources optimally from 
Mobility Management Entity (MME) to base stations (BS)s 
sectors. Finally, we present the simulation results for the per¬ 
formance of our rate allocation algorithm. 

Index Terms —Optimal Resource Allocation, Frequency Reuse, 
Utility Proportional Fairness 

I. Introduction 

In recent years, the cellular system is witnessing a massive 
increase in real-time applications. This increase needs to 
be modeled mathematically for providing efficient resource 
allocation and scheduling methods. The problem of resource 
allocation for cellular systems has been under investigation 
for the last two decade. However, most of the research work 
conducted focuses on concave utility functions for achieving 
global optimal solutions. Otherwise, the achieved solution is 
not optimal. The concave utility functions lack to describe 
real-time applications as shown in m. But the sigmoid utility 
function provides a very good representation of these appli¬ 
cations. However, most of the work conducted on sigmoid 
utility functions provide sub-optimal or approximate solutions 
for resource allocation optimization problem, e.g. El, El 
provides approximate allocation solution and g) curve fits the 
sigmoid functions to the closest concave functions to solve 
with conventional methods. Others only solves for sigmoid 
utility function above a certain threshold to ensure concavity 
of the function 131, |6l. 

In addition, most of the prior work on proportional fairness 
resource allocation optimization considers proportional fair¬ 
ness over rates 0, 0. Nowadays, there is a big variety in the 
applications running under cellular systems operators therefore 
fairness should be application-based not rate-based. The best 
representation of application-based fairness is sigmoid utility 


for real-time applications and logarithmic utilities for delay- 
tolerant applications. In ||8l- lfTol . the authors present utility 
proportional fairness for a cellular system with a single carrier. 
The authors use logarithmic and sigmoid utility functions 
to represent delay-tolerant and real-time applications, respec¬ 
tively. In 0, the rate allocation algorithm gives priority to 
real-time applications over delay-tolerant applications when 
allocating resources as the utility proportional fairness rate 
allocation policy is used. 

In this paper, we include both sigmoid and concave utility 
functions for a cellular system with multiple cells. We show 
how to formulate the resource allocation optimization problem 
with frequency reuse into a convex optimization problem. 

A. Our Contributions 

Our contributions in this paper are summarized as: 

• We introduce a novel rate allocation optimization problem 
with frequency reuse that solves for utility functions that 
are logarithmic and sigmoid. 

• We show that the proposed optimization problem is con¬ 
vex and therefore the global optimal solution is tractable. 
In addition, we present a distributed rate allocation algo¬ 
rithm to solve it. 

The remainder of this paper is organized as follows. Section 
HU presents the problem formulation. Section [III] shows opti¬ 
mality and duality of proposed problem. In Section [IV] we 
present our distributed rate allocation algorithm with frequency 
reuse for the utility proportional fairness optimization problem. 
Section[V]discusses simulation setup and provides quantitative 
results along with discussion. Section [VI] concludes the paper. 

II. Problem Formulation 

We consider a cellular network model that include cells with 
sectors. Therefore, we include frequency reuse. In Figure Q] 
we consider LTE mobile system consisting of K BSs in K 
cells each cell is divided into L sector (e.g. 3 sectors) and M 
UEs distributed in these cells. The rate allocated by the I th 
sector BS to i th UE is given by r\ where i = {1,2 
l = {1,2, ...,2)}. Each UE has its own utility function Ui{ri) 
that corresponds to the type of traffic being handled by the i th 
UE. Our objective is to determine the optimal rates that the I th 
carrier BS should allocate to the UEs. The utility functions are 
given by Ui{n) = Ui{rj + rf + ... -(- rf ) where l r \ = r i 



Fig. 1. System Model. 


and r = {ri,r 2 , ...,rM}- The optimization problem for utility 
proportional fairness: 

M 

max Ui(rl + r? + ... + r\) 

i=l 

L M 

subject to ^ r\=ri , ^r] < R l , ... 

1=1 i= 1 * ' 

M L 

J2^<R L , Y,R l = R, 

i= 1 Z=1 

r{>0, J = 1,2,...,L,*=1,2,...,M. 


where f? z is the allocated rate by the MME to the I th sector, 
and R is the sum of the allocated rates to all sectors. We 
assume that the same frequency band is allocated to all sectors 
therefore avoiding interference. So we have the assumption 
that R l is fixed for all BSs. We assume that UE can’t exist 
in two sectors simultaneously, i.e. if r\ ^ 0 then r\ = 0 for 
q^l. 

We assume the utility functions Ui(r\ + rf +... + rf) to be 
a concave or a sigmoid functions. In our model, we use the 
sigmoid utility function, as in 0, fTTl . given by 


Ui(r\ + rf 


0=Ci( 


1 _|_ e -Oi(Ef=l r \~ b i) 


-di (2) 


where Cj = ' and di = ■ We use the logarithmic 

utility function, as in DU, fl3l . given by 


Ui{r\ + rf + ...+rf) 


log(l + Eti r D 

log(l -T k,r Tnax ) 


(3) 


where r max is the 100% rate and ki is the rate of increase.. 


III. Optimality and Dual Problem 
In the optimization problem (Q}, since the objective func¬ 
tion argmax Ui(r] + rf + ... + rf) is equivalent to 

arg max Jfili log(Ui(rf+rf + ...+rf)), and given that log of 


sigmoid function is concave as shown in j8j then optimization 
problem © in convex and there exists a tractable global 
optimal solution. 

The primal problem in © can be converted into the dual 
problem following similar steps as in (8). The UE and BS 
separate optimization problems follow from similar steps as 
in f8|, fl4l . The utility proportional fairness in the objective 
function of the optimization problem © is achieved by the 
MME which ensure fair allocation to the sectors to achieve 
equal shadow price for all users. This is done by setting pi = p 
for all l. We define the aggregated bids of the I th sectors by 
W l = w i( n ) an so we have W l = pR l and by summing 
over l we have 


J2w l =pJ2n l =pR 


1=1 


1=1 


Then it follows from W l = pR l and @ that 

R l (n) = — - R 

Ef=i 

which is computed in the MME and guarantees fairness. 


(4) 


(5) 


IV. Distributed Optimization Algorithm 

The distributed resource allocation algorithm is an imple¬ 
mentation of UE and BS dual problems, see i), M for 
more details, and ©. Our algorithm allocates resources from 
multiple cells simultaneously with utility proportional fairness 
policy. The algorithm is divided into the i th UE algorithm 
shown in Algorithm © the I th sector algorithm shown in 
Algorithm ©, and MME algorithm shown in Algorithm ©. 
In Algorithm ©, © and ©. the i th UE starts with an 
initial bid wj(l) which is transmitted to the I th sector. The 
I th sector sends the aggregated bids from all UEs under its 
coverage W l (n) = Yfi=i w i( n ) to MME. MME calculates 
the difference between the received aggregated bid W l (n) 
and the previously received bid W (n — 1) and exits if it is 
less than a pre-specified threshold 6 for all sectors. Otherwise, 
MME calculates the sector rates R l (n ) = .E. R and 

sends it to the corresponding sectors. The I th sector uses 
R l {n) to calculate the shadow price pi(n) = — i=1 I fl’ i< ' n ' > and 
sends its value to all the UEs in its coverage area. The i th 
UE receives the shadow prices pi from all sectors that covers 
it and calculates the new bid. This process is repeated until 
| W l (n) — W l (n — 1)| is less than the threshold S. 

V. Simulation Results 

Algorithm © , © and © were applied to various loga¬ 
rithmic and sigmoid utility functions with different parameters 
in MATLAB. The simulation results showed convergence to 
the global optimal solution. In this section, we present the 
simulation results of the users placed in Figure |2] The utility 
function parameters for each user is shown in Table U In the 
following simulations, we set 6 = 10~ 3 . 
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Algorithm 1 The i th UE Algorithm 
Send initial bid wf 1) to BS 

loop 

Receive shadow price p(n) from BS 
if STOP from BS then 

Calculate allocated rate r° pt = 

else 

Calculate new bid Wi(n) = p(n)ri(n) 
Send new bid Wi(n) to BS 

end if 
end loop 


Algorithm 2 The I th sector Algorithm 

Receive sector rate R * l ( 0) from MME {Let R l ( 0) = 

loop 

Receive bids w\{n) from UEs {Let w\(ff) = 0 Vt} 
if | w\(n) — w\{n — 1)| < 8 Vi then 
Allocate rates, H° pt = Wi j n ) to i th UE 

STOP 

else 

Calculate aggregated bids W l (n) = Y^iLi w l( n ) 
Send new aggregated bids W l (n) to MME 
Receive sector rate R l (n) from MME 
Calculate pi(n) = ^0 
Send new shadow price pi (n) to all UEs 
end if 
end loop 


A. Rates for 50 < R < 1150 

In Figure [2 we show the optimal rates of users in the 

1 st sector versus BS rate R. The optimal resource allocation 
is content-aware. The users with real-time application (i.e. 
sigmoid utilities) are allocated resources first. In real-time ap¬ 
plications allocation, the user with the steepest utility function 
(largest a) is allocated first as shown in Figure [3] 

B. Bids for 50 < R < 1150 

In Figure [4] we show the optimal bids of users in the 1 st 
sector versus BS rate R. The users’ bids per resource increases 




Fig. 3. The rates allocated r \ from sector 1 in the three cells. 


Algorithm 3 MME Algorithm 

loop 

Receive bids W l (n) from I th sector {Let W l ( 0) = 0 Vi} 

if \W l (ri) — W l [n — 1)| <8 Vi then 

Allocate sectors rates, R l ° pt = , to I th sector 

Er=i W l (n) 

STOP 

else 

Calculate R\n) = 

V ’ Ei = i W l (n) 

Send new sectors rates R l (n) to I th sector 

end if 
end loop 


as the available resources for allocation in sectors are more 
scarce, i.e. small values of R , as the number of users is fixed. 
Therefore, we have a traffic-dependent pricing. Provided this 
traffic-dependent pricing, the network providers can flatten the 
traffic specially during peak hours by setting traffic-dependent 
bandwidth resource price, which gives an incentive for users 
to use the network during less traffic hours. 

VI. Conclusion 

In this paper, we formulated a novel rate allocation opti¬ 
mization problem with frequency reuse. We considered mobile 
users running two different types of applications, i.e. real¬ 
time and delay-tolerant applications, with utility proportional 
fairness allocation policy in cellular networks. We proved that 
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Fig. 4. The bids wj from sector 1 in the three cells. 


TABLE I 

Users and their utilities 


Sector 1 BS A 

A1 

Sig a = 3, b = 10.0 

A4 

Log k = 1.1, Tmax — 100 

A2 

Sig a = 3, b = 10.3 

A5 

Log k = 1.2, Tmax — 100 

A3 

Sig a = 1, b = 10.6 

A6 

Log k = 1.3, Vmax — 100 

Sector 2 BS A 

A7 

Sig a = 3, b = 10 

A10 

Log k = 1, Tmax — 100 

A8 

Sig a = 3, b = 11 

All 

Log k = 2, rmax — 100 

A9 

Sig a = 1, b = 12 

A12 

Log k = 3, r-max — 100 

Sector 3 BS A 

A13 

Sig a = 3, b = 15.1 

A16 

Log k = 10, Tmax : 100 

A14 

Sig a = 3, b = 15.3 

A17 

Log k = 11, Tmax : 100 

A15 

Sig a = 3, b = 15.5 

A18 

Log k = 12, Vmax — 100 

Sector 1 BS B 

B1 

Sig a = 3, b = 10.9 

B4 

Log k — 1.4, Tmax — 100 

B2 

Sig a = 3, b = 11.2 

B5 

Log k = 1.5, rmax — 100 

B3 

Sig a = 1, b = 11.5 

B6 

Log k = 1.6, Tmax = 100 

Sector 2 BS B 

B7 

Sig a = 3, b = 13 

B10 

Log k — 4, Tmax — 100 

B8 

Sig a = 3, b = 14 

Bll 

Log k = 5, Tmax = 100 

B9 

Sig a = 1, b = 15 

B12 

Log k = 6, r-max = 100 

Sector 3 BS B 

B13 

Sig a = 3, b = 15.7 

B16 

Log k — 13, Vmax -— 100 

B14 

Sig a = 3, b = 15.9 

B17 

Log k = 14, rmax = 100 

B15 

Sig a = 3, b = 17.3 

B18 

Log k = 15, Tmax — 100 

Sector 1 BS C 

Cl 

Sig a = 3, b = 11.8 

C4 

Log k — 1.7, Tmax — 100 

C2 

Sig a = 3, b = 12.1 

C5 

Log k = 1.8, Tmax — 100 

C3 

Sig a = 1. b = 12.4 

C6 

Log k = 1.9, rmax = 100 

Sector 2 BS C 

Cl 

Sig a = 3, b = 16 

CIO 

Log k = 7, rmax = 100 

C8 

Sig a = 3, b = 17 

Cll 

Log k = 8, rmax = 100 

C9 

Sig a=l, b = 18 

C12 

Log k = 9, rmax = 100 

Sector 3 BS C 

C13 

Sig a = 3, b = 17.5 

C16 

Log k = 16, rmax = 100 

C14 

Sig a = 3, b = 17.7 

C17 

Log k = 17, rmax = 100 

C15 

Sig a = 3, b = 17.9 

C18 

Log k = 18, Tmax — 100 


the global optimal solution exists and is tractable for mobile 
stations running delay-tolerant and real-time applications. We 
presented a distributed algorithm, running on UE, BS and 
MME, for allocating resources to mobile users in different 
cells and sectors optimally. Our algorithm ensures fairness in 
utility percentage achieved by the allocated resources for all 
users. Therefore, the algorithm gives priority to the users with 
adaptive real-time applications with providing a minimum QoS 
for all service subscribers. We showed through simulations that 
our algorithm converges to the optimal rates. 
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